The atomic-start description of NiO* 
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We have calculated magnetic properties and the electronic structure of NiO both in the paramag- 
netic and in magnetically-ordered state as well as zero-temperature properties and thermodynamics 
within the strongly-correlated crystal-field approach. It is in agreement with a Mott's suggestion 
that NiO is an insulator due to strong electron correlations. We have quantified crystal-field, spin- 
orbit and magnetic interactions of the Ni 2+ ion in NiO. We have obtained that E^d ^ Eof(=2.0 
cV)>E s „ o (=0.29 eV)»E mQ9 (=0.07 eV). The orbital moment of 0.54 \i B amounts at K, in the 
magnetically-ordered state, to about 20% of the total moment (2.53 (Mb)- Our studies indicate that 
it is the highest time to "unquench" orbital magnetic moment in 3d solid-state physics and the 
necessity to take always into account strong intra-atomic correlations among d electrons and the 
intra-atomic spin-orbit coupling. 

PACS numbers: 75.25.+Z, 75.10.Dg 
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NiO attracts large attention of the magnetic commu- 
nity by more than 50 years. Despite of its simplicity (two 
atoms, NaCl structure, well-defined antiferromagnetism 
(AF) with Tn of 525 K) and enormous theoretical and 
experimental works the consistent description of its prop- 
erties, reconciling its insulating state with the unfilled 3d 
band is still not reached 
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The aim of this paper is to report the calculations of 
the magnetic moment, specific heat and the low-energy 
electronic structure of bulk NiO within the conventional 
ionic picture completed with inter-site spin-dependent 
interactions. We attribute this moment and this low- 
energy electronic structure to the Ni 2+ ions in the rhom- 
bohedrally distorted NaCl structure. The approach used 
can be called the quasi-atomic approach as the starting 
point for the description of a solid is the consideration 
of the atomic-like electronic structure of the constitut- 
ing atoms/ions, fn the present case we start from the 
atomic-like electronic structure of eight electrons in the 
3d shell of the Ni 2+ ions in the oxygen octahedron. A 
macroscopic body of NiO with the NaCl structure is built 
up from such face-sharing octahedra. 

We have treated the 8 outer electrons of the Ni 2+ ion as 
forming the strongly-correlated atomic-like electron sys- 
tem 3d 8 . Its ground term is described by two Hund's 
rules yielding 5=1 and L=3, i.e. the ground term 3 F [ll|. 
In a solid such the localized strongly-correlated electron 
system interacts with the charge and spin surroundings. 
The charge surrounding has predominantly the octahe- 
dral symmetry owing to the NaCl-type of structure of 
NiO. Our Hamiltonian for NiO consists of two terms: 



the single-ion-like term Hd of the 3d 8 system and the 
Hd-d intcrsitc spin-dependent term. Calculations some- 
how resemble those performed for rare-earth systems, 
see e.g. Ref. 0] an d they have been recently applied 
successfully to 3d compounds 0, . For the calcula- 
tions of electronic states of the quasi-atomic single-ion- 
like Hamiltonian of the 3c? 8 system we take into account 
the crystal-field interactions of the octahedral symmetry 
and the spin-orbit coupling (the octahedral CEF param- 
eter E>4=+21 K (=+1.81 meV), the spin-orbit coupling 
constant A s _ = -480 K (= -41 meV) [Tl| ) - It is very 
important to know, that the positive sign of comes 
directly from ab initio calculations for the oxygen anion 
octahedral surroundings. We have chosen the value of B4 
of 21 K in order to reproduce experimentally observed an 
absorption peak at 1.06-1.13 eV 14], sec also Fig. 5 of 
Ref. [8| . The single- ion states under the octahedral crys- 
tal field and the spin-orbit coupling have been calculated 
by consideration of the Hamiltonian |15| : 



Hd = B±(0 { 1 
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The charge-formed fine electronic structure, shown in 
Fig. 1, contains three groups of localized states. The 
higher groups are at 1.06-1.14 eV and 1.91-2.05 eV. For 
low- and room-temperature properties the lowest triplet, 
originating from the cubic subterm 3 A2 g , is the most 
important as the higher states are not thermally pop- 
ulated. The existence of this triplet in the spin-orbital 
space indicates that the Ni 2+ ion is unstable to lattice off- 
octahedral distortions, that splits the triplet. The value 
of B4= +21 K corresponds to another characteristics of 
the cubic crystal field Dq of 1260 K (=875 cm" 1 ). His- 
torically for a configuration with a F ground term (d 2 , 
d 3 , d 7 , d 8 ) the energy separation between 3 A 2s and 3 T2 g 
subterms is denoted as lODq. It should be noted that 
in theoretical approaches different values for lODq are 
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FIG. 1: The calculated fine electronic structure of the highly- 
correlated 3d 8 electronic system in the paramagnetic state, a) 
the 21-fold degenerated 3 F term given by two Hund's rules: 
5=1 and L=3. b) the effect of the cubic octahedral crystal- 
field with Bi=+21 K. c) the combined action of the octahe- 
dral crystal field and the spin-orbit coupling with A s _ = -480 
K. d) the effect of a trigonal distortion - the splitting of the 
ground triplet amounts to 12 K for E>2= +50K. The positive 
B° yields the singlet lower - such a situation is favored by the 
Jahn- Teller theorem. 



considered like only 0.5 eV in Ref. 0, For getting 
the proper insulating gap an enormous splitting between 
spin-up and spin-down e g states of 10.8 eV was intro- 
duced in Ref. |4|, that according to us is not physically 
justified. 

For the description of the trigonal distortion observed 
in NiO the cubic Hamiltonian (1) is necessary to trans- 
form in order to have the z quantization axis along the 
cube diagonal. In case of the cubic NaCl structure the 
trigonal distortion is realized by stretching or compress- 
ing along the cube diagonal. The CEF Hamiltonian to 
describe such the distortion takes the form: 



H d = -2/35 4 (C>2 - 20V2Ol) + B 2 0° 2 



(2) 



The obtained eigenfunctions of the lowest spin-orbital 
triplet, for 5 4 =+21 K and A s _ = -480 K yielding 
A s _ o /10Dq=0.038, take a form (only \L Z ,S Z ) components 
of the full \L, S, L Z ,S Z ) function are shown): 

V>i = 0.740 |0,+1) -0.522|+3,+l) + 0.418|-3,+l) 

ip = 0.744 |0,0) -0.467|+3,0) + 0.467 |-3,0) 

= 0.740 |0,-1) +0.522|-3,+l) - 0.418|+3,+l> 

For these functions ^=±0.9945 and L z =±0.297, what 
give the total moment of m z =±2.2883 hb- The L z value 
of 0.297 is close to the first-order spin-orbit correction 
(Ref. [HI, p. 450) of 8-A so /10Dq=8-0. 038=0.30. 

For the trigonal distortion for B% >0 the triplet splits 
into lower singlet and higher doublet. For B\ <0 the 
doublet is lower. For £?2=+50 K the splitting of the 
ground triplet amounts to 12 K. 
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FIG. 2: The calculated temperature dependence of the fine 
electronic structure of the highly-correlated 3d 8 electron sys- 
tem in the magnetically-ordered state below Tjv of 525 K. In 
the paramagnetic state the electronic structure is tempera- 
ture independent unless we do not consider, for instance, a 
changing of the CEF parameter due to the thermal lattice 
expansion. 



The magnetic field, external or internal in the case 
of the magnetically-ordered state, polarizes two states 
of the doublet, as is shown in Fig. 2. The intersite 
spin-dependent interactions Hd-d cause the (antiferro- 
)magnetic ordering. They have been considered in the 
mean-field approximation with the molecular-field coeffi- 
cient n acting between magnetic moments rrid=-(L+g e S) 
- 2.002324. The value of n in the Hamiltonian 



H, 



d-d 
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has been adjusted in order to reproduce the 
experimentally-observed Neel temperature. The fit- 
ted value of n has been found to be -11.55 meV/ 
H B (— -200 T/fis)- It means that the Ni ion in the 
magnetically-ordered state experiences the molecular 
field of 503 T (at K). This field as well as the magnetic 
moment become smaller with the increasing tempera- 
ture and vanish above Tjy. The calculated value of the 
magnetic moment at K in the magnetically-ordered 
state amounts to 2.53 /is, Fig. 2. It is built up from the 
spin moment m s of 1.99 \ib (52=0.995) and the orbital 
moment m of 0.54 [is. Fig. 3. The increase of m in 
comparison to the paramagnetic state, ± 0.26 /ib, is 
caused by the further polarization of the ground-state 
eigenfunction by the internal molecular magnetic field. 
The orbital moment is quite substantial being about 
20% of the total moment. The calculated by us moment 
at 300 K amounts to 2.26 /is {m s = 1.78 /is, m = 
0.48 fiB ) nicely reproducing the experimental result of 
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FIG. 3: The calculated temperature dependence of the Ni 2+ - 
ion moment in NiO. At K the total moment mjv; of 2.53 {ib 
is built up from the orbital m D and spin m s moment of 0.54 
and 1.99 fis, respectively. 

2.2±0.3 \i B for the Ni moment at 300 K 16]. It should 
be noted that this novel moment value is remarkably 
larger than old results of 1.64-1.70 /is, which have been 
generally up to now used in theoretical calculations 
[3, Q ■ The magnetic X-ray experiment of Ref. 0] has 
revealed a substantial orbital moment of 0.32 /j,b and 
the spin moment of 1.90 [Ib at 300 K. 

The trigonal distortion is important for the detailed 
formation of the AF structure and the direction of the 
magnetic moment but it only slightly influences the spin 
and orbital momenta values. The trigonal distortion 
yielding the singlet charge-formed ground state leads to 
the moment direction lying in the plane perpendicular to 
the cube diagonal. Exactly such the moment direction is 
observed in NiO. Actually, the magnetic ordering occurs 
along the (112) direction within this diagonal perpendic- 
ular plane due to a further slight distortion within the 
(111) plane @. 

The present model allows, apart of the ordered mo- 
ment and its spin and orbital components, to calculate 
many physically important properties like temperature 
dependence of the magnetic susceptibility, temperature 
dependence of the heat capacity (shown in Fig. 4), the 
spectroscopic g factor, the fine electronic structure in 
the energy window below 3 eV with at least 20 localized 
states, Fig. 1. The spike-like peak in c(T ) at 7jv is in 
perfect agreement with experimental data |17| obtained 
on a single-crystal specimen that yields " very large, very 
narrow peak of 65 cal/Kmol" [lg- Ref- U3 provides a 
detailed critical analysis of literature experimental data 
with the tabulated recommended data, that are shown in 
Fig. 4. As the best description of the c(T) dependence 
in the range 0<T<250 K they give the Debye value of 8 
of 580 K. In Fig. 4 we provide our rough summation of 
the calculated c<j(T) contribution and the lattice contri- 
bution with a value of 8 of 650 K. We find the overall de- 
scription of c(T) as remarkably good owing too so large 
temperature interval and the A- type peak at T^r. The 
magnetically-ordered state of NiO has lower energy than 
the paramagnetic one by 3.4 kJ/mol (= 35 meV/ion) at 
K, Fig. 2. Although the molecular field acting on the Ni 
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FIG. 4: The calculated temperature dependence of the 3d 
contribution Cd(T) to the heat capacity of NiO. The spike 
peaks to 214 J/Kmol. A value of 6 R, R is a gas constant, for 
the high temperature Dulong-Petit heat is shown. Points are 
the recommended data by Ref. 18] after a critical analysis of 
literature experimental data. 



moment is so large as 503 T it splits CEF states, however, 
without their reorganization in the lowest triplet. The 
splitting of the spin-up and spin-down states amounts to 
1560 K (135 meV). We think that this energy excitation 
among the lowest triplet has been revealed by Newman 
and Chrenko almost 50 years ago [l4|, who has found 
a 240 meV excitation. From temperature dependence 
they have concluded that it is connected with the anti- 
ferromagnetic ordering, exactly as comes out from our 
approach. In our picture this excitation decreases with 
increasing temperature vanishing at Tjv, see Fig. 2. 

We would like to point out that our approach should 
not be considered as the treatment of an isolated ion 
only - we consider the Ni 2+ ion in the oxygen octahedron. 
The physical relevance of our calculations to macroscopic 
NiO is obvious - the NaCl structure is built up from the 
face sharing Ni 2+ octahedra. In the perfect structure 
all Ni ions experience the same charge (= crystal field) 
and spin interactions, what means that all Ni ions have 
the same electronic structure. The good reproduction 
of macroscopic properties proves that all atoms equally 
contribute to these properties. But one should remember 
that due to a finite size of each sample always some atoms 
on the surface, say 1-2 %, will have another electronic 
structure like it was discussed in Ref. 

In contrary to a general conviction within the mag- 
netic community about the adequateness of the strong 
crystal-limit to 3d-ion compounds we have found a re- 
markably good description of NiO in a rather weakly- 
intermediate CEF limit. The octahedral crystal field is 
25 times stronger than the spin-orbit coupling but it does 
not break intra-atomic correlations among 3d electrons 
(after giving up two electrons during the formation of 
a solid NiO, from the 4s states). Such a physical sit- 
uation is the basis for a developed Quantum Atomistic 
Solid State Theory (QUASST) for 3d-atom containing 
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compounds [Igj. We treat our QUASST approach as a 
continuation of studies of Van Vleck on correlation of 
macroscopic properties with discrete crystal-field states 
of transition-metal ions [2(j • 

In conclusion, we have calculated magnetic properties 
and the electronic structure of NiO both in the paramag- 
netic and in magnetically-ordered state as well as zero- 
temperature properties and thermodynamics. We have 
quantified crystal-field (the leading parameter B4 = +21 
K), spin-orbit (-480 K, i.e. like in the free ion [Tj) and 
magnetic interactions (B moJ of 503 T and n= -200 T/[i B ). 
In our approach E dd > E c _f(=2.0 eV)>E s _ o (=0.29 
eV)»E maff (=0.07 eV). The orbital and spin moment 
of the Ni 2+ ion in NiO has been calculated within the 
quasi-atomic approach. The orbital moment of 0.54 \xb 
amounts at K in the magnetically-ordered state, to 
about 20% of the total moment (2.53 hb)- Despite of 
using the full atomic orbital quantum number L=3 and 
5=1, the calculated effective moment from the tempera- 
ture dependence of the susceptibility amounts to 3.5-3.8 
/iB, i.e. only 20 % larger value than a spin-only value 



of 2.83 hb- Our success is related to the fact, that we 
take into account the existence of very strong correla- 
tions among electrons. Good description of many phys- 
ical properties indicates that 3d electrons in NiO are in 
the extremely strongly-correlated limit, i.e. in the ionic 
limit confirming a priori our atomic-start assumption. In 
our approach we can trace the Jahn- Teller effect, in the 
spin-orbital space, and the breaking of the time reversal 
symmetry, equivalent to the appearance of the ordered 
magnetic state, at the atomic scale. 

We can mention that within the same approach we 
have described a few 3d compounds like FeBr2 [l3|. 
LaCo0 3 , Na 2 V 3 7 and recently CoO HJ. The local- 
ized crystal-field states discussed in the present paper 
are already in a technical use in a ultrafast manipulation 
of the antiferromagnetism of NiO • 

* dedicated to John Van Vleck and Hans Bethe, pio- 
neers of the crystal-field theory, to the 75 th anniversary 
of the crystal-field theory, and to the Pope John Paul II, 
a man of freedom in life and in Science. 
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